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PREFACE 


Recent  studies  in  the  application  of  lifting  surfaces  and  pointed 
slender  cones  to  both  glide  vehicles  and  advanced  re-entry  vehicles 
have  focused  attention  upon  the  aerodynamic  performance  of  such  con¬ 
figurations  in  sustained  hypersonic  flight  through  the  lower  ionosphere. 
In  addition,  interest  has  been  exhibited  by  aerodynamicists  and  heat- 
transfer  engineers  in  the  fluid-mechanical  and  thermal  phenomena  which 
occur  in  the  leading-edge  region  of  a  sharp  wing  which  is  in  hypersonic 
flight  at  even  lower  altitudes. 

This  Memorandum  examines  a  theoretical  problem  which  is  highly 
relevant  to  both  of  these  applications:  the  analysis  of  the  low-density 
hypersonic  flow  near  the  leading  edge  of  a  sharp  flat  plate.  The  anal¬ 
ysis  should  be  applicable  to  the  determination  of  the  incipient  effects 
of  intermolecular  collisions  upon  the  flow  of  a  rarefied  gas  in  the 


transition  regime. 
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SUMMA.RY 


The  initial  interactions  between  a  hypersonic  free  stream  and 
molecules  emitted  (reflected  diffusely)  by  a  cold  flat  plate  immediately 
downstream  of  its  sharp  leading  edge  are  analyzed  by  near-free-molecule 
theory.  The  results  give  the  variations  of  the  pressure,  shear,  and  heat 
transfer  at  the  wall.  These  are  found  to  increase  downstream  from  the 
leading  edge  and  reach  a  maximum  plateau,  and  a  subsequent  decay  of 
the  wall-transfer  properties  to  continuum  boundary- layer  values  can 
be  inferred.  The  magnitude  and  location  of  this  plateau  compare  well 
with  available  viscous  pressure  interaction  data.  An  indication  of  the 
initial  growth  rate  of  the  disturbed  merged  layer  is  also  obtained 
from  the  analysis,  it  is  found  that  this  region  is  wedgelike  and  has 
a  slope  of  approximately  45  deg  for  all  hypersonic  Mach  numbers  and 
for  wall-to-free-s tream- temperature  ratios  of  the  order  of  unity. 
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SYMBOLS 

numerical  proportionality  constants 

proportionality  factors  relating  molecular  model 
characteristics  to  flow  parameters 

skin-friction  coefficient  and  Stanton  number 

exponential  integral 

empirical  constant  (essentially  the  inverse  thermal- 
accomodation  coefficient) 

velocity-distribution  function 

collisfon-density  distribution  (Eq.  (7))  of  molecular 
group  i 

wall-to-free-stream  temperature  factor  (Eq.  (22)) 

Bessel  function  of  the  second  kind 

mean  free  path  of  a  molecular  group  i  diffusing  into  the 
atmosphere 

free  path  of  a  molecular  group  i  diffusing  into 
molecules  j  only 

ratio  of  mean  free-stream  velocity  relative  to  the 
plate  to  the  average  (directed)  thermal  velocity  in 
the  free  stream  (Eq.  (24)) 

Mach  number,  ratio  of  mean  relative  velocity  to  the  speed 
of  sound  in  the  free  stream 

flux  of  molecules  per  unit  area  (molecules/sec,  unit  area) 
surface  pressure 

Reynolds  number  based  on  free-stream  conditions 

mutual  collision  cross  section 

temperature  (deg  abs) 

magnitude  of  the  vector  velocity 

cross-stream  and  streamwise  velocity  components 
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=  free-stream  mean  (directed)  thermal  agitation 

=  streamwise  and  cross -stream  coordinates  (origin  at 
leading  edge) 

=  nondimensional  coordinates  x/L  and  y/L 

e  J  e 

=  normalization  factor  (Eq.  (5)) 

=  scattering  angles  for  primary  encounters 
=  collision  rate  per  unit  volume 

=  polar  coordinates  of  a  point  (0  measured  relative  to 
the  normal  to  the  plate) 

=  normal  momentum  coefficient 

=  Maxwellian  mean  free  path  in  the  free  stream 
=  molecular  number  density  (molecules/unit  volume) 

=  functions  describing  wall  emission 
=  hypersonic  interaction  parameter  V Re 

SUBSCRIPTS 

=  molecules  emitted  (reflected)  from  the  surface 

=  molecules  resulting  from  collisions  of  molecule  e 
with  n,,  and  so  on 

=  free-molecule  (collisionless)  flow 
=  generic  subscripts  denoting  classes  of  molecules 
=  incident  molecules 
=  free-stream  molecules 
=  surface  of  the  plate 
=  free-stream  stagnation  conditions 
=  free-stream  static  conditions 

=  indicates  averaging  over  the  distribution  of  x 


-1- 


I.  INTRODUCTION 
(12  3) 

Recent  measurements  ’  ’  '  of  the  surface  pressure  on  flat  plates 
in  hypersonic  flow  at  low  densities  have  shown  that  there  exists  a 
measurable  length  downstream  of  the  leading  edge  over  which  the 
concepts  of. first-order  boundary- layer  interaction  theory  do  not 
apply.  There  is  evidence  of  a  pressure  plateau,  the  magnitude  of 
which  depends  both  on  the  Mach  number  and  also,  apparently,  on  the 
temperature  of  the  flow  relative  to  the  wall  temperature.  Downstream 
of  this  plateau  the  pressure  decays  at  a  rate  which  is  in  qualitative 
agreement'  with  boundary- layer  interaction  concepts  but  which  does 
not  agree  with  currently  available  first-order  theories.  Upstream' 
of  it  the  pressure  must. rise,  a  phenomenon  which  is  entirely  outside 
of  the  scope  of  the  continuum- flow  analytical  model.  Schlieren 
photographs^^  show  this  initial  growth  of  the  disturbed  region  and 
suggest  an  exponential  shape  for  it.  It.  is  interesting  to  note  that  . 
this  structure  of  the  flow  is  independent ' of  the  density,  provided 
only  that  the  leading  edge  is  sufficiently  sharp  to  prevent  the  form¬ 
ation  of  a  bow  wave.  The  flow  passes  through  all  . regimes  of  aero¬ 
dynamics,-  from  f ree-molecule  to  continuum,  in  a  certain  distance 
downstream  from  the  leading  edge. 

In  this  Memorandum,  the  origin  of  the  flow  over  a  semi- infinite 
plate  which  has  a ' sharp  leading  edge  aligned  with  the. stream  is  in¬ 
vestigated  by  a  kinetic- interaction  technique.  The  model  is  illus¬ 
trated  schematically  in  Fig.  1.  The  incident . molecules,  which  belong 
to  a  molecular  group  defined  by  the  characteristics  of  their  mean 


Fig.  I  —  Sketch  of  initial  interactions 
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motion,  approach  the  plate  nearly  tangentially.  (The  analysis  is 
restricted  to  high  Mach  numbers.)  Another  group,  particles  emitted 
by  the  wall,  diffuse  into  this  atmosphere  at  right  angles  to  the 
free  stream.  These  particles  collide;  this  collision  knocks  them 
out  of  their  groups  and  replaces  them  by  particles  of  a  new  group, 
the  characteristics  of  which  can  be  derived  from  the  dynamics  of  the 
encounter.  The  paths  of  these  first- scattered  molecules  are  shown 
on  Fig.  1  by  the  darker  lines;  the  figure,  and  the  analysis  as  well, 
is  restricted  to  the  case  in  which  the  mean  velocity  of  the  emitted 
particles  is  low  compared  to  that  of  the  incident  ones,  relative  to 
a  coordinate  system  fixed  to  the  plate. 

The  spatial  distribution  of  the  first  collisions  (zone  of  first 
collisions)  is  calculated.  It  is  characterized  by  the  convective 
distance  between  collisions,  a  relative  mean  free  path  of  one  group 
diffusing  into  the  atmosphere  of  the  other  group--not  by  the  mean 
free  path  of  either  group  by  itself.  Obviously,  if  the  mean  velocities 
of  each  group  are  different,  the  zone  of  first  collisions  will  be 
geometrically  distorted. 

By  knowing  thus  the  source  distribution  of  the  first- scattered 
molecules,  their  paths  can  be  calculated  and  their  decay  due  to  second 
collisions  or  their  capture  by  the  surface  can  be  found.  The  scale 
of  the  zone  of  second  collisions  is  characterized  again  by  a  relative 
mean  free  path  of  the  scattered  molecules  in  the  atmosphere  of  the 
others;  this  mean  free  path  is  different  from  that  of  the  first 
collisions.  An  examination  of  the  total  normal  momentum  due  to 
molecules  of  the  various  groups  captured  by  the  surface  yields  the 
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local  surface  pressure,  which  is  the  primary  result  of  this  analysis. 

The  analysis  of  this  interdiffusion  and  scattering  of  successive 
groups  of  molecules,  assuming  a  one-sided  emitting- collecting  surface, 
can  by  the  present  technique  be  continued  past  the  first- scattered 
group  as  long  as  distinct  characteristic  interaction  parameters  for 
each  group  (mean  free  paths)  can  be  defined.  This  turns  out  to  be 
practical  when  there  are  only  a  few  families  in  a  given  region  in 
space  and  when  the  differences  in  their  characteristics  are  large. 

In  other  words,  the  analysis  is  practical  in  the  zone  of  initial 
collisions  (up  to  the  third)  when  the  flow  is  hypersonic  (nearly 
Newtonian)  and  when  the  wall  temperature  (the  mean  velocity  of  the 
emission  by  the  wall)  is  low  compared  to  the  free- stream  stagnation 
temperature.  Nonetheless,  significant  results  are  obtained  because 
the  very  first  interaction  between  a  free- stream  molecule  and  one 
emitted  by  a  surface  brings  about  a  violent  change  in  the  contribution 
of  the  former  molecule  to  the  normal  momentum  transfer  to  the  surface: 
The  collision  scatters  them  at  an  average  of  45  deg  to  the  relative 
direction  of  impact.  Thus  it  transforms  one-half  the  free-stream 
momentum  per  molecule  into  normal  momentum  directed  at  the  surface 
(cold  surface) . 

Downstream  of  the  initial  collision  zone  one  distinguishes  a 
transition  zone  in  which  further  collisions  bring  about  the  final 
establishment  of  an  equilibrium  shear  layer.  In  an  equilibrium  shear 
layer,  the  transverse  momentum  transfer  from  the  free  stream  to  the 
wall  proceeds  by  a  collision-diffusion  process  rather  than  by  direct 
molecular  convection.  The  present  analysis  cannot  be  carried  out 


formally  into  the  transition  zone,  although  some  of  its  character¬ 
istics,  based  on  the  indications  of  the  first- collision  calculation, 
are  presented  in  the  text.  Figure  1  shows  schematically  the  structure 
of  the  leading-edge  zone  to  the  formation  of  a  continuum  boundary 
layer  as  it  is  visualized  in  this  discussion. 

Similar  initial-collision  calculations  performed  for  blunt-body 

(4  5) 

flows'  ’  have  shown  that  in  spite  of  their  restrictions,  the  theory 
yields  a  physically  meaningful  picture  of  the  phenomena  and  leads  to 
the  definition  of  important  interaction  criteria  and  similarity  para¬ 
meters.  This  is  the  principal  purpose  of  the  present  analysis.  We 
view  it  as  a  subtle  form  of  dimensional  analysis  rather  than  as  a 
quantitative  theory.  We  hope  to  obtain  from  it  a  qualitative  frame¬ 
work  useful,  perhaps,  for  constructing  a  semiempirical  correlation  of 
the  data,  which  at  present  escapes  generalization.  (Note  the  wide 
scatter  of  the  experimental  measurements  when  plotted  against  the  con¬ 
tinuum  interaction  parameter.)  Accordingly,  we  make  several  simplify¬ 
ing  assumptions  in  the  calculations,  provided  that  they  do  not  distort 
the  essential  physical  phenomenon.  An  accurate  solution  of  the  leading 

edge  problem  will  have  to  be  attempted  along  different  lines,  perhaps 

/£\ 

by  use  of  a  moment  method  of  solving  the  Boltzmann  equations.  An 
important  feature  of  the  analysis  will  be  its  ability  to  span  the  en¬ 
tire  transition  from  free-molecule  to  continuum  flow. 

It  is  perhaps  worth  mentioning  in  this  discussion  the  difference 
between  the  present  technique  and  other  methods  which  have  also  been 
called  first-collision  analyses — for  instance,  the  method  proposed 
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(7  8) 

by  Willis.  ’  They  are  not  equivalent.  Willis  finds  a  first-order 
correction  to  the  free- stream  distribution  function  by  iterating  on  a 
simplified  (Krooks 1  model)  Boltzmann  integral,  starting  with  the 
same  two  Maxwellian  groups  of  molecules,  free- stream  molecules  and 
those  emitted  by  the  wall.  This  technique  differs  from  the  present 
method  in  that  the  mean- free- path  characteristic  of  the  interaction 
is  not  a  priori  a  constant.  The  present  iteration  (not  in  a  mathe¬ 
matical  sense)  on  successive  groups  of  scattered  molecules,  although 
numerically  more  approximate,  avoids  the  need  for  computer  techniques 
and  yields  a  clearer  physical  picture  of  the  field  of  the  interactions 
than  does  Willis'  method. 
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In  order  to  describe  the  interdiffusion  of  the  molecular  groups 
we  consider  the  length  of  flight  of  a  particle  of  the  group  i  in  an 
atmosphere  of  particles  j 


ij  ©.  .  P  .S.  .1  V.  +  V. 

ij  3  lj'-i  -j 


where  0. .  is  the  collision  rate  of  i  with  j.  S. .  is  the  mutual  collision 

ij 

cross  section,  and  |  V.  +  V.j  is  the  magnitude  of  the  relative  (vector) 

i  1 

precollision  velocity. 

In  general,  the  factors  in  Eq.  (1)  are  functions  of  both  velocity 
and  physical  spaces.  A  mean  free  path  can  be  obtained,  in  principle, 
by  carrying  out  the  appropriate  integrations.  In  particular,  the 
local  decay  of  beam  i  by  collisions  with  j  is 


dn.  = 

i 


P.S.  .  V.  +  V. 

1  ii'-i  -v 


n.  f  ,e 

j.  i 


dVdedC 


V  9  C 


where  f.  is  the  distribution  function  of  molecules  i  emitted  from  a 

l 

point  in  the  polar  direction  9  at  velocity  V^,  and  where  the  integra¬ 
tion  is  performed  over  all  velocity  and  physical  spaces  (polar  co¬ 
ordinates  0,  C)  . 

Equation  (2)  leads  to  a  difficult  implicit  problem  unless  the 
product  p  |  +  Vj |  is  constant  or  prespecified  in  geometric  space. 

If  j  is  to  represent  a  particular  group  (for  instance,  the  free- stream 
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molecules),  then  such  an  assumption  restricts  the  validity  of  the 

calculation  to  the  extreme  near- free-molecule  zone  and  cannot  correctly 

•k 

reproduce  the  interesting  feature  of  the  leading- edge  problem.  How- 
ever,  we  can  consider  j  as  representing  all  the  groups  in  the  zone  of 
collisions  with  the  groups  i.  This  assumption  requires  that  the 
magnitude  of  the  relative  precollision  velocity  vector  be  adequately 
represented  by  an  average  of  all  the  possible  collisions.  Now,  we 
shall  concern  ourselves  only  with  the  initial  collisions  which  follow 
an  encounter  between  a  molecule  emitted  (reflected)  from  the  wall 
(type  e)  and  a  free-stream  molecule  (type  n) .  We  shall  also  assume 
V£  <<  V  ,  that  is,  that  the  wall  emission  is  cold  in  comparison  with 
the  stagnation  temperature  of  the  flow.  In  this  case  the  relative 
vector  velocity  for  the  most  probable  spherically  elastic  first 
collision  is  simply  V  ,  for  all  possible  second  collisions  (1 /-v/TjV^, 
and  for  nearly  all  third  collisions  1/2  (see  Fig.  2).  Only  about 
7  per  cent  of  the  quaternary  molecules  begin  to  have  an  upstream 
velocity  component. 

It  follows,  provided  that  there  exist  relatively  well-defined 
space  domains  populated  only  by  molecular  groups  which  have  collided 
not  more  than  three  times,  that 

P.V.+V.  =  P.lv.+V.l  =  Const  for  the  group  (3) 

j  -1  -j  jl-x  -jl 

is  physically  significant.  It  turns  out  that  such  domains  exist  as 

'  Vs 

This  assumption  is  implied  in  previous  analyses  of  blunt-body 
flow.J^’5) 

See,  for  instance,  Kennard ,  Kinetic  Theory  of  Gases,  Chapter  3 
(McGraw) .  The  vector  diagrams  are  based  on  momentum  conservation  for 
a  two-dimensional  gas. 
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Fig.  2 — Most  probable  scatterings  in  successive  collisions 
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long  as  the  molecular  motion  is  .predominately  unidirectional  in  the 
initial- collision  zone  at  the  leading  edge  for  hypersonic  flow  over 
a  cold  surface. 

The  foregoing  is  a  fundamental  restriction  on  the  validity  of 
the  present  near- free-molecule  technique.  The  other  assumptions 
which  are  made  in  this  calculation  and  which  are  outlined  below  are 
not  essential  and  can  be  removed  if  desired.  Their  only  purpose  is 
to  permit  a  closed- form  solution. 

The  individual  particle  velocity  V.  of  a  group  i  is  replaced  by 
an  average  ^  over  velocity  space.  Note  that  the  dependence  of 
on  velocity  always  feeds  back  through  the  intermediary  of  particle- 
particle  collisions  to  the  description  of  the  distribution  in  velocity 
space  of  the  fundamental  group  of  molecules  e  emitted  by  the  wall 
in  a  given  polar  direction  6,  about  which  nothing  is  known.  Converse¬ 
ly,  the  distribution  of  the  average  velocity  in  the  polar  directions 
9  from  the  collision  center  can  be  conserved  without  added  difficulty 
by  writing 


Y 


(©) 


For  example,  in  the  case  of  Lambert  diffuse  emission  from  a 
Maxwellian  oven  into  a  half- space,  we  have 


Tt  TT 

¥(0)  =  cos  e,  -  -  <  e  <+  2 


as  usual.  Since  is  not  a  function  of  G,  the  beam- decay  formula, 
Eq.  (2),.  reduces  to  the  simple  form 
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Vc@) 


V(e) 


dC 


(4) 


which  upon  integration  through  the  polar  space  ©  from  the  source  re¬ 
duces  to  the  classical  form  applicable  to  a  homogeneous  molecular 
* 

field;  is  here  a  constant  mean  length  characterizing  the  collisions 

of  that  group  with  all  other  molecules . 

We  shall  also  assume  in  the  calculations  that  the  motion  of  all 
molecules  is  restricted  to  two  dimensions  in  a  plane  parallel  to  the 
mean  flow  and  perpendicular  to  the  leading  edge.  It  is  clear  that 
in  the  application  to  a  two-dimensional  problem  (infinite  span)  this 
assumption  does  not  change  the  pattern  of  the  interactions  at  all 
and  affects  the  numerical  constants  only  slightly. 


THE  EMISSION  FROM  THE  SURFACE 

Let  be  the  polar  density  distribution  of  the  ng  molecules 

per  unit  surface  and  time  emitted  by  an  element  dx^  of  the  plate  at 
Xp  into  a  wedge  +  A  6  about  ©.  By  using  Eq.  (4),  the  density  of 
their  collisions  in  space  (collisions  per  unit  time,  per  element  of 


volume  2A©£dQ  is 


c 

n 

e 


1  c®)  ^ 


(5) 


where  is  a  normalization  factor 


*A  refinement  worth  noting  is  this:  The  proper  value  of  V.  which 
enters  Eq.  (4)  in  the  exponent  is  not  numerically  identical  to  £he 
"mean"  velocity  in  the  sense  of  temperature  of  that  molecular  group. 
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e  +  Ae 


l  e  -  Ae 

2Ae 


*«»  de 


* 


In  order  to  obtain  the  total  number  of  collisions  per  unit  volume 
and  time  at  a  generic  point  in  the  half- space  above  the.  plate  it  is 
necessary  to  integrate  along  the  emitting  surface.  Thus 


0)  e 


(e)  V 


where  the  coordinates  C  and  6  (9  measured  from  the  perpendicular  to 
the  plate)  relate  a  point  in  space  to  the  source  of  emission,  and 
where  n^  =  ne^xj  is  generally  a  function  of  position  along  the  sur¬ 
face  as  dictated  by  a  mass  balance  among  the  emitted  and  all  the 
intercepted  molecules  at  x^ .  It  may  include  a  mass  source,  such  as 
sublimation  of  the  plate  material.  Clearly,  this  depends  implicitly 
on  the  results  of  the  analysis,  and  it  is  necessary  to  accept  either 
an  approximate  or  an  iterative  procedure  to  continue.  A  relatively 
simple  scheme  is  to  expand  ng  in  series  in  x^  with  coefficients  to 
be  determined  later  from  auxiliary  mass- conservation  conditions. 
However,  it  was  found  that  although  this  is  conceptually  simple,,  it 
complicates  the  subsequent  integrals  sufficiently  to  demand  a  numerical 
solution. 

The  functions  0  (which  is  equivalent  to  1  )  and  Y  describe  the 
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geometric  character  of  wall  emission.  Since  very  little  is  known  about 
this  character,  we  must  make  some  assumptions  in  order  to  facilitate 
the  integration.  Some  typical  models  considered  are  shown  in  Fig.  3. 
These  models  include  the  classical  Lambert  Law  (Model  I),  a  density 
distribution  of  a  specular  type  based  on  indications  of  the  recent 
experimental^^  and  theoretical  studies,  and  the  simplest  possible 
model  of  a  beam- like  monochromatic  emission  from  the  wall,  of  an  ar¬ 
bitrary  degree  of  specularity  described  by  0  (Model  V)  . 

Equations  (8),  (9),  and  (10)  below  show  the  6- integrals  (Eq.  (7)) 
for  three  of  the  models.  The  result  for  Model  I  is  approximate  and 

_x^  1 

is  obtained  by  replacing  the  exponential  e  by  a  factor  - r— 

(1  +  xZ) 

which  behaves  similarly  to  it.  Using  only  a  constant  term  in  the  ex¬ 
it 

pansion  for  n^ 


Model  V 


G  L 


e  e 


n 

e 


e 


P 

cos  0 


(10) 


Figure  4  shows  lines  of  constant  parameter  (G^L^/n^)  in  the 

*K  is  the  Bessel  function  of  the  second  kind  and  Ei  is  the  ex¬ 
ponential  integral. 
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f  =  COS  0 
©  =  cos  0 


LAMBERT  LAW 


+  =  cos  0 

©  =  SIN  20 


COSINE  VELOCITY,  SPECULAR 


*  =  I 

©  =  COS  0 


UNIFORM  VELOCITY,  DIFFUSE 


*  =  I 

©  =  SIN  20 


UNIFORM  VELOCITY,  SPECULAR 


3 L 


MONOCHROMATIC  BEAM  AT  0 


Fig.  3  —  Wall- emission  models 
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initial -interaction  zone  for  Models  I  and  II.  Model  V  fails  to  re¬ 
produce  the  two-dimensional  nature  of  the  leading- edge  zone  and  the 
upstream  propagation  of  the  disturbance  seen  in  Fig.  4.  However,  the 

extent  of  this  zone  is  of  the  order  L  ,  and  it  will  be  shown  later 

e’ 

that  the  length  characterizing  the  molecular  exchanges  with  the  surface 
(induced  pressure)  is  of  the  order  L^M.  Therefore,  when  the  Mach 
number  is  sufficiently  high,  the  stand-off  distance  of  the  interaction 
zone  relative  to  the  leading  edge  and  the  influence  of  the  type  of 
the  wall  emission  become  unimportant. 

Note  that  Eqs.  (8)  and  (9)  already  Imply  the  assumption  of  a 
constant  (streamwise)  surface- emission  density.  However,  the  Gg 
function  for  Model  V,  Eq.  (10),  permits  this  x- dependence  to  be  carried 
on.  This  may  offer  advantages  in  further  studies,  such  as  the  inclu¬ 
sion  of  mass  addition,  which  outweigh  the  shortcomings  of  this  simple 
distribution. 

THE  DISTRIBUTION  OF  FREE- STREAM  MOLECULES 

The  space  density  of  the  incident  molecules  diminishes  as  they 
penetrate  the  interaction  zone  because  of  collisions  which  knock  them 
out  of  their  class.  Two  molecules  of  another  class  per  collision, 
the  characteristics  of  which  can  easily  be  obtained  from  the  collision 
dynamics,  are  generated  in  their  place.  The  rate  of  disappearance  of 
incident  molecules,  n,  due  to  collisions  with  the  wall-emitted  ones, 
e,  is,  per  unit  volume  and  unit  time 

G 

e 


P 

n 

P 


(ID 
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which  is  also  twice  the  number  of  en,  the  first- scattered  molecules, 
generated  per  unit  volume  and  time  at  that  point.  Since  the  motion 
of  en  is  known,  it  is  conceptually  possible  to  calculate  Ggn  in  space 
and  the  rate  of  disappearance  of  n  by  collisions  with  en,  and  so  on. 
The  local  density  of  the  n  molecules  in  space  is  given  by  the  solution 
of  the  conservation  equation 


dp 

Mn 

un 


+  v 


n  dy 


.1 


.  G  . 
J  nj 


(12) 


In  this  general  form  the  problem  is  extremely  complicated.  We 
shall  consider  in  what  follows  only  the  first  iteration;  that  is, 
we  shall  disregard  all  but  the  first  collision.  This  can  be  inter¬ 
preted  physically  by  stating  that  we  are  not  distinguishing  between 
the  two  tertiary  molecules  (for  example  enn  and  ene)  and  the  two 
molecules  which  generated  them.  From  the  standpoint  of  the  cross¬ 
stream-momentum  flux,  which  is  the  principal  result  of  this  analysis, 
this  is  correct  because  downward-directed  momentum  is  generally 
conserved  in  such  collisions.  From  the  standpoint  of  the  space 
distribution  of  these  molecules,  the  approximation  is  acceptable  be¬ 
cause  the  most  probable  vector  paths  of  the  enn  or  the  ene  molecules 
do  not  differ  too  much  from  those  of  the  generating  ones;  also,  their 
speeds  and  their  mean  free  path  in  the  gas  are  not  greatly  different. 

If  the  flight  Mach  number  is  sufficiently  high,  the  cross- stream- 

* 

diffusion  term  (pseudo-Newtonian  flow)  can  be  neglected.  Then  the 


*A  better  but  more  complicated  solution  can  also  be  obtained  if 
the^secon^  term  in  Eq.  (12)  is  replaced  by  a  self-diffusion  term 
X  a  Pn/dy  .  It  is  interesting  to  note  that  the  result  then  depends 
also  on  the  free- stream  mean  free  path  X. 
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local  density  of  the  incident  molecules  in  the  interaction  zone  is 
simply 

p  =  P  e  (13) 

n(x,y) 


The  restriction  to  hypersonic  flows  is  consistent  with  previously 
discussed  assumptions  concerning  the  existence  of  a  characteristic 
interaction  distance,  L  .  The  density,  P,  in  the  exponent  of  Eq.  (13) 
is  the  total  local  density  (all  groups),  and  the  product  PV^  is  con¬ 
stant  (see  Eq.  (3))  (Va  «  V^) . 

The  exponent 

form  of  sums  of  integral  functions  for  all  the  models  considered. 

The  results  are,  however,  too  complicated  for  further  closed- form 
analysis  except  in  the  case  of  Model  V.  In  that  case,  from  Eq.  (10) 

.  B 

/  \  nr 

P 


"  ^  e 

-  /< 


Ggdx  in  Eq.  (13)  can  be  obtained  in  the 


cos  0 


In 


p3  =  '  (“  ~  P  tan  6  )  (  pv“ 


(14) 


Here  9  is  the  angle  of  the  beam  inclination. 


SPACE  DISTRIBUTION  OF  FIRST -SCATTERED  MOLECULES  (MODEL  V) 

Each  (pn/p  G  )  collision  per  unit  volume  and  time  produces  two 
en  molecules  at  that  point  in  space.  Their  motion  relative  to  the 
collision  center  can  be  specified  in  terms  of  the  motions  of  the 
colliding  pair  without  difficulty  but  at  the  cost  of  considerable 
algebraic  complexity.  The  problem  can  be  simplified  without  essen¬ 
tial  physical  error  by  considering  only  the  most  probable  collision, 
in  which  case  the  two  en  molecules  are  emitted  in  monochromatic 
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beams  from  the  collision  center  along  two  distinct  angles,  7^  and 

7„ •  In  particular,  in  the  extreme  case  when  V  is  zero  (stationary 
^  £ 

scatterer),  7^  -  7^  =  rr/4,  the  vector  velocity  of  both  en  molecules 
is  (l//2)VnJ 

and  the  velocity  components  normal  and  tangential  to 

the  flow  of  each  molecule  are  (1/2)V  . 

n 

However,  a  better  relation  for  the  properties  of  the  scattered 
molecules,  which  maintains  at  least  approximately  their  dependence 
on  the  emission  velocity,  is 


(15) 

(16) 


The  positive  and  negative  signs  refer  to  the  upward  and  the  downward 
family,  respectively. 

The  concern  with  this  new  group  of  molecules  brings  forth  a  new 
characteristic  mean  free  path.  From  the  basic  definition 


L 

en 


V 

en 


V 

— en 


+  V. 


\ 


(17) 


where  the  second  equality  is  obtained  from  the  most-probable- collision 
diagram  and  Eq.  (16).  It  is  notable  that  the  mean  free  path  of  scattered 
molecules  equals  the  Maxwellian  mean  free  path  in  the  free  stream. 

It  is  identical  for  both  the  en-o  and  the  en-e  encounters  (j  in  Eq . 

(17)  is  n  and  e,  respectively)  in  this  particular  case  of  the  cold  wall. 
Thus  this  single  characteristic  length  adequately  describes  the 
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interactions  in  the  domain  of  the  first- scattered  molecules. 

Let  us  consider  separately  the  two  scattered  beams.  We  consider 
first  the  upward  flux  in  order  to  outline  approximately  the  character 
of  the  interaction  zone. 


where  the  first  factor  in  the  integral  represents  the  local  source 
density  of  en  molecules  at  a  generic  point,  and  the  exponentials 

describe  their  decay  along  their  path,  which  is  assumed  to  be  char¬ 
acterized  adequately  by  a  single  mean  free  path,'"L  .  The  integration 
is  performed  from  an  adequate  lower  limit  (the  plate  surface),  denoted 
by  C .  The  source  term  is 


(19) 


where  G£  represents  all  the  collisions  of  molecules,  e,  at  the  point 

C  at  which  the  fraction  P  /p  generates  the  en  molecules  under  con- 
o  n 

sideration. 

For  the  case  9=0  and  Model  V  one  finds  for  a  >  g 


where  K  is  related  to  the  wall-to-free-stream  temperature  ratio 


and  M'  is  a  molecular  speed  ratio  defined  by  Eq.  (24).  This  integral 
is  approximated  somewhat  by  neglecting  an  exponential  factor  in  the 
integrand.  Because  of  this,  it  is  not  valid  in  the  neighborhood  of 
o:  =  P.  An  expression  similar  to  Eq.  (21),  and  valid  for  or  <  P,  is 
also  found  but  for  the  sake  of  brevity  will  not  be  shown. 

Lines  of  constant  density  ratios  of  both  the  n  and  the  en  mole¬ 
cular  groups  (for  two  specific  cases)  are  plotted  in  Fig.  3.  This 
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can  be  compared  with  the  schematic  history  of  the  molecular  motion 
outlined  in  Fig.  1 , 

Downward  Flux 

The  second  family  of  the  en  molecules  is  directed  toward  the 
surface.  Noting  the  relative  magnitude  of  L^,  which  measures  the 
distance  from  the  surface  of  the  plate  from  which  they  originate, 
and  L  ,  which  measures  the  length  of  their  path,  we  conclude  that 
in  the  cold-wall  case  the  overwhelming  majority  of  these  en  molecules 
will  strike  the  surface.  Moreover,  the  most  probable  secondary  col¬ 
lisions  (en-e  or  en-o)  do  not  change  the  normal  momentum  intercepted 
by  the  surface.  Only  tertiary  collisions  ene-o,  and  so  on,  transform 
further  a  portion  of  the  streamwise  momentum  into  momentum  normal  to 
the  plate.  If  the  calculation  is  to  be  interrupted  by  this  family  of 
secondary  molecules,  then  a  better  result  for  the  total  momentum  inter¬ 
change  with  the  surface  is  obtained  by  neglecting  the  tertiary  collisions 
entirely  and  assuming  that  all  the  en^  molecules  generated  in  the 
interaction  zone  reach  the  plate.  The  corresponding  integral,  analogous 
to  Eq.  (18),  does  not  contain  the  exponential- decay  factor.  Carrying 
out  the  integration  with  the  appropriate  limits  (over  the  quadrant 
above  the  plate)  one  finds  to  the  same  approximation  the  following 
expression  for  the  flux  (molecules/second-unit  length)  to  the  plate 


This  is  shown  in  Fig.  4  for  various  Mach  numbers  (and  7^  =  T1/4,  as  before). 
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DISCUSSION 


The  analysis  has  resulted  in  several  characteristic  parameters 

of  interest.  Since  the  number  of  scattered  molecules ,  n  ,  must  be 

1  e 

equal  to  the  flux  of  molecules  to  the  plate,  and  since  in  this 
approximation  n^  was  assumed  to  be  constant  and  therefore  equal  to 
the  free- stream  cross -flow  flux 


n 

e 


Pv„ 


n  n 


(24) 


where  vT  is  the  average  directional  free- stream  thermal  agitation. 
Thus  the  first  parameter  which  describes  the  density  of  the  free¬ 
st  ream  group  involves  a  factor  proportional  to  the  reciprocal  of  the 
Mach  number. 

Using  this  in  Eq.  (14),  it  is  seen  that  the  decay  of  the  free- 

_  Oi 

stream  molecules  occurs  in  a  rectangular  region  (P  /M'),  which 

is  equivalent  to  (L^  •  ML^) .  It  is  stretched  in  the  flight  direction 
by  M.  As  M  goes  to  infinity  (Newtonian  flow),  the  whole  interaction 
zone  (the  shock  layer  and  the  boundary  layer)  lies  identically  on 
the  surface. 

The  well-known  relations  between  the  interaction  length,  L^, 
and  gasdynamic  parameters  is 


where  X  is  the  mean  free  path  in  the  free  atmosphere,  and  T  is  the 
wall  temperature.  It  follows  that  the  physical  distance  from  the 
leading  edge,  x,  in  which  a  fixed  percentage  (or  all)  the  free- stream 
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molecules  have  collided  is  not  a  function  of  Mach  number  but  of  gas 
density  and  of  the  wall-to-free-stream  temperature  ratio  T  /T^  only. 

There  is  obviously  a  singularity  for  Tw  =  0  which  has  no  physical 
meaning. 

The  field  of  decay  of  the  en  group  is  more  complex  than  that  of 
the  n  group.  It  depends  on  the  spatial  density  of  e-n  collisions 
and  on  the  decay  characteristics  of  the  en  group  itself.  The  decay 
occurs  in  a.  streamwise  length  of  the  order  of  I^M,  which  is  the  dimension 
of  the  first- collision  zone  generating  the  en  molecules,  and  in  a 
cross-stream  distance  zone  which  has  thus  the  dimensions 

T 

X  -2  •  X 

It  is  approximately  square  and  on  the  order  of  the  free- stream  mean 
free  path.  It  is  only  slightly  stretched  streamwise  with  the  square 
root  of  the  wall- to- free- stream  temperature  ratio,  and  it  is  essentially 
independent  of  the  Mach  number  (see  also  Fig.  5).  This  appears  to 
be  a  fundamental  characteristic  of  the  initial  growth  of  the  disturbed 
layer  over  the  plate.  It  has  a  slope  of  the  order  of  unity  rather 
than  the  slope  of  infinity  which  is  implied  by  continuum  boundary- layer 
theory.  It  does  not  seem  to  follow  the  exponential  function  suggested 
in  Ref.  1 . 

The  flux  of  secondary  molecules  to  the  plate  increases  from  the 
leading  edge  to  a  maximum,  as  seen  in  Fig.  6.  Beyond  this  point  it 
decreases  because  the  interaction  zone  is  being  depleted  of  the  free- 
stream  molecules,  n,  which  generate  the  flux.  It  is  of  course  incorrect 
to  imply  that  the  total  mass  and  momentum  flux  to  the  surface  decreases; 
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Calculated  local  density  of  incident  and 
scattered  molecules 
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tertiary  and  further  molecules  of  all  kinds  make  up  the  difference. 

We  shall  assume  that  the  initial-collision  zone  is  bounded  on  the 
downstream  side  by  the  locus  of  maxima  of  the  flux  of  en  molecules. 

In  this  paper  we  wish  to  interrupt  the  "iteration"  with  the  family 
of  en  molecules,  yet  we  should  like  to  sketch,  at  least  approximately, 
the  further  development  of  the  flow.  Following  the  initial-collision 
zone  is  a  transition  zone  in  which  the  molecular  families  multiply 
in  a  geometric  progression,  and  there  are  no  definite  space  domains 
in  which  one  family  or  a  reasonably  small  number  of  families  dominates. 
This  can  be  seen  from  consideration  of  the  simple  vector  diagrams  of 
the  successive  collisions.  Already  one  of  the  tertiary  collisions 
(en-e)  generates  a  molecule  which  moves  upstream.  Thus  the  direction¬ 
ality  of  the  motion  which  makes  this  kinetic  near- free-molecule  analysis 
feasible  is  rapidly  lost  in  the  transition  region,  and  more  complete 
integral  kinetic  methods  must  be  used  to  study  it. 

The  normal  momentum  exchange  with  the  surface,  which  jumps  from 
zero  to  one-half  the  free- stream  momentum  per  molecule  upon  the  initial 
collision,  reaches  a  maximum  at  the  end  of  this  zone;  it  varies  little 
in  the  transition  region,  because  the  further  conversion  of  stream- 
wise  momentum  into  cross- stream-directed  momentum  is  roughly  a 
geometric  progression  in  terms  of  the  number  of  collisions.  Gross 
calculations  for  the  density  of  tertiary  collisions  were  carried  out, 
and  these  indicated  that  indeed  the  total  mass  and  the  corresponding 
momentum  flux  to  the  plate  extend  smoothly  (asymptotically)  past 
the  maximum  of  the  en  flux  into  the  transition  region. 


It  is  clear  that  the  normal  momentum  flux  to  the  surface  must 
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reach  a  plateau.  The  foregoing  arguments  indicate  that  this  is  of 

the  order  of  the  flux  existing  at  the  edge  of  the  initial  collision 
* 

zone.  We  shall  assume  this  in  the  following  discussion  of  the 
induced  pressure. 


Note  that  it  was  assumed  that  all  the  (en)^  molecules  reach  the 
wall.  It  was  shown  that  tertiary  collisions  are  thereby  approxi¬ 
mately  accounted  for  in  what  concerns  the  downward  momentum  flux 
intercepted  by  the  surface. 
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III.  THE  INDUCED  PRESSURE 


A  number  of  physically  significant  flow  properties  (wall  shear, 
velocity  profiles,  second  approximation  to  the  local  density,  etc.) 
can  be  derived  from  the  previous  calculations.  In  this  Memorandum 
we  shall  discuss  primarily  the  problem  of  the  pressure  at  the  surface. 
This  is  the  only  quantity  that  has  been  measured  under  conditions 
which  fit  approximately  the  assumptions  of  this  analysis  and,  therefore, 
the  one  for  which  some  verification  of  the  results  can  be  provided. 

The.  most  recent  and  the  most  complete  series  of  tests  is  that  re¬ 
ported  by  Nagamatsu  and  Sheer. ^  The  measurements  are  reproduced  in 
Fig.  7.  The  major  portion  of  the  measurements  are  well  into  the 
continuum  regime;  however,  the  curves  do  show  a  well-defined  plateau 
near  the  leading  edge.  According  to  what  was  stated  previously, 
this  maximum  occurs  in  the  transition  regime,  and  its  magnitude  equals 
approximately  the  pressure  at  the  edge  of  the  initial- collision  domain. 

To  calculate  the  wall  pressure  we  must  consider  the  free-stream 
molecules  (Eq.  (14))  contributing  m  v^,  v^,  to  the  normal  momentum 
per  unit  surface  plus  the  en  molecules  (Eqs.  (16)  and  (23)),  carrying 

approximately  1/2  n  mV  (V.  -  V  /V  )  momentum  per  unit  surface, 
r  J  en  e  1  e  n 

and  finally  the  reaction  due  to  the  emission  from  the  wall,  ng  m  Vg. 

The  expression  for  the  induced  pressure  becomes 


We  shall  begin  by  studying  the  induced -pressure  maximum  in  order 
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B.  COLD  TESTS  TS=2300°R;  Tw<v540°R 


FIG.  7  —  MEASURMENTS  OF  INDUCED  PRESSURE01 
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to  adjust  the  numerical  constants  (which  were  approximated  throughout 
the  analysis)  and  to  verify  the  truth  of  the  model  we  formulated. 

We  remark  that  the  pressure  maximum  is  expected  to  occur  at  a,  value 
of  a  larger  than,  but  on  the  order  of,  unity.  However,  a/M'  is 
expected  to  remain  small  because  the  Mach  number  must  be  large. 
Accordingly,  expanding  the  exponential  in  Q'/m'  and  carrying  the 
(analysis  out  to  first  order  in  that  parameter,  we  find  that  the  max¬ 
imum  of  Eq.  (26)  occurs  at  o  given  by  the  solution  of 

max 


_ 2 

M'2  1  -  V  /V 

e  n 
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max  -  1  1 

M'  ,) 
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Therefore 
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where  the  successive  forms  of  Eq.  (28)  represent  successive  numerical 
approximations  which  are  allowable  for  a/M'  small  and  M'  large. 

By  the  same  degree  of  approximation,  Eq.  (26)  can  be  expanded  to  read 
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Combining  Eqs.  (28)  and  (29),  the  expression  for  the  maximum  induced 
* 

pressure  becomes 


By  comparing  this  approximate  form  to  the  full  Eq.  (26),  we  find  an 

error  of  less  than  5  per  cent  at  most  in  the  range  of  the  experimental 

data  (M  >  8)  . 
o 

DISCUSSION  AND  COMPARISON  WITH  DATA 

Before  proceeding  in  a  comparison  of  the  results  of  this  theory 

to  the  available  experimental  points,  it  is  necessary  to  introduce 

some  proportionality  constants,  omitted  in  the  analysis,  by  the 

assumption  of  a  two-dimensional,  nearly  Newtonian  molecular  model. 

The  nonlinear  form  of  Eq.  (30)  in  both  independent  variables  M'  and 

V  /V  makes  this  difficult  to  do  by  purely  graphical  means.  However, 
e  n 

a  physical  analysis  of  the..  mo-del...pxQvidfiS-  the  .fallowing,  guidex. — 


*The  free-stream  pressure,  P^,,  is  defined  herein  in  terms  of  the 
momentum  exchange  with  an  element  of  surface  parallel  to  the  flow  in 
equilibrium  with  it;  that  is 

Poo  *  (p  vT)  +  ne  vT  «  2Pvtvt> 
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In  the  course  of  the  analysis  three  separate  molecular  families 
were  defined,  and  their  interactions  were  described  in  terms  of  the 
characteristic  velocities  of  the  families.  Each  of  these  characteris¬ 
tic  velocities  is  an  unspecified  average  over  the  velocity-distribution 
function  of  the  family,  and,  moreover,  in  each  case  the  motion  is 
assumed  to  be  two-dimensional  (neglecting  the  component  parallel  to 
the  leading  edge)  . 

We  are,  therefore,  led  to  introduce  the  following  proportionality 
cons tants : 

a.  The  free-stream  family  is  described  by  the  parameter  M',  which 
is  a  sort  of  molecular  speed  ratio  based  on  the  mean  directed  thermal 
velocity  of  the  gas.  We  set  it  proportional  to  the  true  Mach  number 

as  follows 

V 

M'  «  —  «  CM 
—  o 

VT 

b.  The  surface-emitted  family  is  characterized  by  the  mean 
velocity  V£ .  Eliminating  the  free-stream  velocity,  we  write 


Both  Vg  and  v^,  are  averages  over  distributions  of  the  same  character, 
and  their  ratio  needs  no  additional  constants.  However,  is  not 
generally  known.  Instead,  the  surface  temperature  is  known  and  thus 
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where  F  is  essentially  the  inverse  thermal-accommodation  coefficient. 

It  should  be  recalled  that  the  entire  analysis  (with  Model  V)  has 
made  the  assumption  9  =  0y  which  corresponds  to  an  assumption  of  a 
normal  momentum-accommodation  coefficient  of  unity  (diffuse  reflection). 

It  should  also  be  borne  in  mind  that  the  "cold  wall"  assumption 
implied  here  requires  that  ~  1;  and  therefore  not  only 

T^/Tjj  ~  1,  but  also  F  must  be  on  the  order  of  unity. 

Combining  these  definitions,  the  velocity  ratio,  V  /V  ,  can  also 
be  expressed  in  terms  of  the  wall-to-stagnation  temperature  ratio 


This  formulation  is  convenient  for  comparison  with  data  such  as  those 

in  Ref.  1,  which  consist  of  series  of  tests  at  constant  T  /T  . 

w  s 

However,  the  use  of  the  ratio  T^/T^  as  a  basic  variable  probably 
would  be  preferable  in  future  analytic  as  well  as  experimental 
work. 

c.  The  scattered  family  is  described  by  the  velocity  V  , 
again  a  two-dimensional  characteristic.  Accordingly,  we  set 
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where  D  absorbs  the  factor  1/2  used  thus  far.  There  is,  of  course, 
a  cross-coupling  between  the  scattered  family  and  the  wall-emission 
temperature . 

We  have,  consequently,  four  basic  constants: 

C,  D=  numerical  proportionality  factors  relating 
to  the  molecular  model 

9,  F  =  surface-interaction  parameters,  essentially 
the  normal  momentum  (a)  and  inverse  ther- 
mal-accoimiodation  coefficients  (a) 

The  first  two  stem  from  the  simplifications  of  this  analysis  and 
could  be  calculated  at  the  cost  of  numerical  effort.  The  last  two 
(in  this  or  in  another  form)  are  the  ever-present  kinetic  surface- 
interaction  parameters  which  would  require  the  assumption  of  a 
surface -interaction  model  such  as  Schamberg '  s  ,  and  this 

remains  the  principal  unknown  in  free-molecule  and  near-free- 
molecule  theory. 

In  the  following  we  shall  make  initially  the  common  assumption 
that,  there  is  perfect  surface  accommodation,  i.e., 

9-0  F  =  1 

and  use  the  two  available  empirical  equations  describing  the  induced 
pressure  (two  series  of  tests  each  at  constant  Tw/Tg)  to  evaluate 
C  and  D.  The  procedure  and  the  evidence  that  F  is  not  unity  in 
these  tests,  as  well  as  remarks  concerning  the  effect  of  this  on 
the  present  comparison  between  the  theory  and  the  experiment,  are 
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given  in  the  Appendix. 

Introducing  these  constants  into  Eq.  (30)  and  rearranging,  we 
obtain 


Using  the  procedure  described  in  the  Appendix  and  the  data  of 
Ref.  1  we  find  the  most  reasonable  estimate  of  the  constants  to  be 

C  -  0.33 

D  »  25 


Equation  (31)  is  shown  graphically  in  Fig.  8  in  a  comparison 
with  the  experimental  measurements .  The  correlation  is  remarkably 


There  is  also  available  data  showing  the  plateau  of  hypersonic 
induced  pressure  obtained  by  J.  Aroesty:  Pressure  Distribution  on  Flat 
Plates  at  Mach  4  and  Low  Density  Flow,  Univ.  of  California  Inst.  Eng. 
Res.  Rep.  HE  150-157,  July  1958.  This  is  not  used  for  comparison  with 
this  theory  because  the  tests  were  conducted  on  an  insulated  model 
with  nearly  stagnation  wall  temperature,  and  the  condition  Ve/Vn  <  <  1 
implied  in  the  analysis  was  definitely  not  met.  Qualitatively,  how¬ 
ever,  there  are  no  contradictions  between  Aroesty's  measurements  and 
the  present  results. 


2.8 


COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 
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good,  since  it  covers  a  large  range  of  both  independent  parameters 
of  the  problem,  which  enter  the  equation  in  a  nonobvious  way.  How¬ 
ever,  the  quantitative  results  should  not  be  overemphasized  for 
several  reasons:  The  molecular  model  is  approximate;  the  experi¬ 
mental  measurements  are  not  certain  (it  is  pointed  out  in  the  Appen¬ 
dix  how  an  even  better  fit  can  be  obtained  if  one  assumes,  for 
example,  a  partial  accommodation  to  surface  conditions);  and  finally, 
some  of  the  test  conditions  do  not  comply  very  closely  with  the  con¬ 
ditions  for  validity  of  the  analysis--a  high  Mach  number  and  a  very 
low  wall  temperature. 

Perhaps  the  most  significant  result  is  the  prediction  of  a  de¬ 
crease  in  the  maximum  induced  pressure  with  an  increase  of  wall-to- 
stagnation  temperature  ratio.  This  trend  seems  definitely  to  be 
verified  by  the  experiment.  On  the  other  hand,  it  contradicts  the 
prediction  of  the  merged-layer  analysis  based  on  slip-free  continuum 
concepts .  In  order  to  bring  this  fact  into  focus,  the  data,  the 

present  kinetic  theory,  and  Oguchi 's  wedge-like-flow  solution 
are  compared  in  Fig.  9.  It  appears  that  they  coincide  for  the 
particular  case  of  =  0.08,  the  only  experimental  test  discussed 

by  Oguchi,  but  when  applied  to  the  low-stagnation-temperature  case 
his  results  diverge  significantly  from  both  the  data  and  the  present 
results . 

It  is  also  interesting  to  observe  that  the  present  theory  pre¬ 
dicts  in  the  high-Mach-number  case  a  linear  variation  of  the  induced- 
pressure  ratio  with  Mach  number 


G.  9  — COMPARISON  OF  DATA  AND  LEADING-EDGE 
INTERACTION  THEORIES  FOR  VARYING  Tw/Ts 
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while  Oguchi's  result  gives  a  parabolic  increase  of  this  ratio 


Lim  Ap 


where  G,  a  function  of  T^/T'  ,  lies  between  1/4  and  tr/4.  The  former 
seems  more  reasonable  than  the  latter  on  the  ground  of  physical  argu¬ 
ments  . 

The  initial  pressure  rise  (first-order  in  a)  is  easily  derived 
from  Eq.  (26).  In  terms  of  the  three  independent  parameters  of  the 
problem,  the  Mach  and  Reynolds  numbers  and  the  wall- to- free- stream 
temperature  ratio,  and  where  A  and  B  are  numerical  constants,  it  is 


At  a  fixed  Mach  number  and  increasing  rarefaction  (Re-*0),  the  induced 
pressure  rise  is  dominated  by  the  retroforce  due  to  the  wall  emission 
(the  second  term);  this  is  in  accord  with  free-molecule  theory.  With 
increasing  density,  that  is,  towards  continuum  flow  (Re-»oo),  the  slope 
of  the  initial  pressure  rise  tends  to  infinity.  However,  the  maximum 
induced  pressure  at  =  amax  remains  finite  regardless  of  density,  as 
it  depends  only  on  the  Mach  number  and  the  wall- to- free- stream  temper¬ 


ature  ratio. 
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*Note  that  P. 
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The  difference  — - - —  is  shown  in  Fig.  10  for  specific  values 

TM. 

1 

of  the  independent  parameters.  Note  that  the  relative  magnitude  as 
well  as  the  location  of  the  maximum  of  these  curves  depends  on  the 
choice  of  the  proportionality  constants,  C  and  D,  and  the  surface- 
accommodation  constant,  F.  The  heat  transfer,  in  particular,  is 
sensitive  to  the  value  of  this  constant.  In  order  to  demonstrate 
this  effect,  two  sets  of  curves  are  shown. 

In  principle,  measurements  of  the  heat- transfer  plateau  such  as 

•k 

are  being  currently  obtained  could  provide  the  third  empirical 
equation  from  which  C,  D,  and  also  F  could  be  determined  simultaneously. 
This  is  discussed  in  more  detail  in  the  Appendix  in  connection  with 
the  preliminary  heat-transfer  data  that  are  available. 

SECOND  APPROXIMATION  TO  n 
_  en 

In  the  foregoing  discussion,  the  density  of  the  surface 

emission  was  assumed  to  be  constant  along  the  plate.  A  second 

iteration  on  n  can  be  carried  out  easily  for  the  initial  zone 
e 

to  give  an  estimate  of  the  effect  of  this  assumption.  To  the 
second  order,  the  local  surface  emission  equals  the  local  flux 
of  the  free  stream,  nn>  and  the  first-scattered  molecules,  n^, 
collected  by  the  surface;  that  is 

“k 

Personal  communications  from  Dr.  H.  T.  Nagamatsu,  General 
Electric  Research  Laboratories,  Schenectady,  N.  Y. 
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(2)  _  (1)  ,  (1) 
n  =nv  +n 
e  n  en 

The  result  is,  for  the  initial  zone  for  which  both  tf/M  <<'  1  and 
-o' 

e  ~  C1  -  a) 


The  second  term  represents  the  increase  of  the  flux  due  to  the  cou¬ 
pling  between  the  surface  emission  and  the  density  of  the  first  scat¬ 
terings  . 

Continuing  as  outlined  for  the  first-order  analysis,  the  pressure 
on  the  surface  downstream  of  the  leading  edge  (x  <  L^)  is 


(34) 
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The  second  approximation  gives  essentially  the  same  form  for 

the  initial  pressure  rise,  which  is  dominated  by  the  product  M'Re, 

(2) 

but  the  slope  of  the  second-order  approximation  (the  constant  Av  ') 
is  higher  by  a  factor  of  about  2. 


SURFACE  SHEAR  NEAR  THE  LEADING  EDGE 

By  following  the  procedure  already  established,  we  can  obtain 

the  first  approximation  to  the  local  shear  coefficient  by  adding 

the  contributions  due  to  the  incident  flux  (mo  V  )  to  those 

n  i  n 

due  to  the  scattered  flux  (mo  v  u  ).  The  emitted  molecules 

en  en  en 

do  not  carry  any  tangential  momentum  in  the  case  of  the  simple 
normal  beam-emission  model  (Model  V)  considered  here. 

We  find  without  difficulty 
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T 

S3  . 

1/2  P  V  2 
n  n 
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O'  ,  -o 
M 1  C 


(35) 


This  expression  reduces  to  the  free-molecule  skin  friction 


C 


f 


FM 


2 
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(36) 


when  Leo  tends  to  infinity  for  any  finite  value  of  x.  This  is 
indeed  the  correct  form  for  the  molecular  models  considered. 

As  before,  simplified  expression  for  the  hypersonic  case  in 
the  initial-collision  domain,  and  for  the  extreme  leading-edge 
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We  find  thus  that  the  skin  friction  behaves  like  the  pressure. 
There  is  an  initial  increase  of  local  friction  downstream  of  the 
leading  edge,  a  maximum,  and  finally  a  decrease  in  the  transitional 
zone.  This  is  due  to  the  flux  of  first- scattered  molecules  to  the 
surface,  which  overrides  the  decrease  of  the  flux  of  the  incident  mole¬ 
cules  and  the  decrease  of  momentum  per  molecule  of  the  scattered  mole¬ 
cules  compared  to  the  incident  ones. 

It  is  notable  that  this  result  is  characteristic  of  the  high-Mach- 
number  case.  Although  this  analysis  is  not  strictly  valid  for  low 
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Mach  numbers  because  of  the  basic  assumptions  discussed  previously, 
the  trend  is  nonetheless  evident.  If  one  sets  M  =  0  in  Eq.  (35) 


At  low  M  the  friction  coefficient  decreases  immediately  downstream 

/io\* 

of  the  leading  edge.  This  is  the  same  result  as  that  of  Liu.  ' 

Liu  neglects  the  change  in  the  space-density  distribution  of  molecules 
brought  about  by  the  interactions  and  therefore  does  not  obtain  the 
term  which  increases  at  all.  In  this  sense  his  physical  model  is 
only  valid  for  the  near-stationary  case;  moreover,  in  the  numerical 
approximation  he  sets  M  <  1  in  his  calculation. 

HEAT  TRANSFER,  SLIP,  AND  TEMPERATURE  Jl)MP 

This  model  is  valid  only  for  the  highly  cooled  plate  (V^  <  <  V  ; 

2 

Tw  <  <M  TJ,  and  therefore  a  recovery  temperature  cannot  be  calculated. 
The  variation  of  the  local  heat-transfer  coefficient  in  the  initial- 
collision  zone  for  a  fixed  wall  temperature  (fixed  V^)  is  of  course 
qualitatively  similar  to  the  variation  of  the  induced  pressure.  Thus, 
we  expect  a  rapid  rise  to  a  maximum  followed  by  a  decay  to  the  value 
of  continuum  theory.  Incidentally,  a  constant  Vg  would  require  a  par¬ 
ticular  heat-absorption  schedule  (or  a  particular  variation  of  ther¬ 
mal  accommodation  to  the  plate)  . 

A  slip  velocity  and  a  temperature  jump  at  the  wall  can  be  calcu¬ 
lated  simply  from  the  foregoing.  The  physical  slip  velocity  for  the 
simple  first-order  model  investigated  here  is 


Liu's  function  G  is  approximately  proportional  to  M  so  that  his 
correction  term  (Re/M^)  G  exhibits  the  same  general  dependence  on  the 
flow  parameters. 
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which  can  be  gotten  directly  from  Fig.  6  or  Eq.  (23).  It  is  seen 
that  in  the  domain  of  validity  of  this  analysis  the  slip  velocity 
decays  from  free  stream  to  about  75  per  cent  of  free- stream  velocity 
at  hypersonic  speeds.  The  decrease  of  is  monotonic  from  the  leading 
edge.  Now,  in  terms  of  slip- flow  theory,  the  wall  shear  is  proportion¬ 
al  to  the  slip  velocity,  which  is  taken  to  be  proportional  to  the 
mass  velocity  gradient  at  the  wall.  The  present  analysis  indicates  an 
inconsistency  in  the  slip-flow  assumption  if  it  is  to  be  extended 
into  near- free -molecule  flow,  because  very  near  the  leading  edge  in 
hypersonic  flow  the  shear  is  shown  to  increase  while  Ug  decreases. 
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IV.  EMPIRICAL  EXTENSION  OF  THE  RESULTS 

The  present  analysis  is  valid  in  the  initial- collision  domain 
only,  but  hopefully  it  can  also  be  used  as  a  framework  for  further 
semiempirical  correlations.  There  is  one  critical  point  on  the 
hypersonic  interaction  curves :  namely,  the  continuum  limit  of  the 
constant-pressure  plateau  (marked  on  the  curves  of  Fig.  11  by  a 
vertical  arrow)  which,  if  it  were  predictable,  would  allow  the  induced 
pressure  to  be  drawn  for  all  values  of  the  flow  parameters. 

A  schematic  induced- pres sure  curve  is  drawn  as  an  insert  in 
Fig.  12.  One  distinguishes  the  leading- edge- pres sure  rise  to  a  hyper¬ 
sonic  interaction  parameter,  xg>  corresponding  to  the  edge  of  the  ini¬ 
tial-collision  region.  Following  this  region,  there  is  the  constant- 
pressure  plateau  between  values  Xg  and  X c  of  the  interaction  parameter. 
This  region  was  previously  identified  with  a  transition  in  which  the 
kinetic  motion  in  the  disturbed  layer  readjusts  itself  to  the  equili¬ 
brium  condition-  characteristic  of  a  continuum  boundary  layer.  From 
this  point  on,  the  further  growth  and  the  displacement  effect  of  the 
boundary  layer  are  governed  by  shear  transfer  and  mass  continuity 
as  postulated  in  the  continuum  analysis  of  interaction  theories. 

The  edge  of  the  initial  collision  zone,  Eq.  (28),  is  determined 
physically  by  the  convective  distance  to  undergo  a  fixed  number  of 
collisions.  The  parameter  which  describes  this  is  the  slip  parameter^^ 
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Fig,  1 1  —  Induced-pressure  measurements  versus 
continuum -interaction  parameter(,) 
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NOTE  : 

TYPICAL  INDUCED 
PRESSURE  VARIATION 


Fig. 12— The  extent  of  the  pressure  plateau 
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where  the  second  equality  is  obtained  from  Eq.  (28).  The  Mach  number 
2 

dependence  ofXg/M  reflects  the  streamwise  distortion  (stretching) 
of  the  zone  of  first  collisions. 

We  postulate  that  the  extent  of  the  transition  zone  is  also 
determined  by  a  fixed  number  of  collisions,  which  are  required  to 
establish  a  new  equilibrium.  Since  the  molecular  motion  in  the  trans¬ 
ition  region  is  predominately  along  45-deg  lines  and  both  the  normal 
and  streamwise  components  of  the  molecular  paths  are  of  the  same  order, 
we  do  not  expect  a  streamwise  distortion  of  this  zone.  The  convective 
distance  required  for  the  fixed  number  of  collisions  is  then  a  con¬ 
stant,  and  we  seek  a  correlation  in  terms  of 


X 


s 


=  Constant 


(42) 


Figure  12  shows  the  results.  (The  use  of  the  Mach  number  as  the 
ordinate  is  purely  a  means  of  displaying  the  data.)  Indeed,  within 
the  accuracy  with  which  the  continuum  edge  of  the  transitional  plateau 
can  be  defined,  Eq.  (42)  is  well  substantiated  with  a  proportionality 
constant  of  about  1.6.  Note  that  the  data  include  a  considerable  range 
of  Mach  numbers  and  an  independent  variation  of  stagnation  temperatures. 

One  is  also  tempted  to  use  the  results  of  this  analysis  to  plot 
velocity  or  density  profiles  through  the  disturbed  layer.  This  can, 
of  course,  be  done  numerically.  However,  we  believe  that  more  pre¬ 
cise  results  than  those  which  can  be  derived  by  an  inspection  of  Fig. 

5  would  not  be  significant;  inasmuch  as  good  arguments  can  be  given 
for  the  validity  of  the  wall- flux  calculation,  the  details  of  tire 
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field  are  probably  grossly  oversimplified  by  this  model.  The  one 
result  worth  noting  is  the  wedge-like  growth  of  the  disturbed  zone; 
a.  sort  of  oblique  characteristic  line  defined  by  a  maximum  in  the 
gradient  of  the  density  of  disturbed  molecules  can  be  seen  to  be 
along  the  line  x  =  y  (hyperthermal  flow),  regardless  of  density. 

This  line  should  not  be  identified,  however,  with  the  continuum  gas- 
dynamic  oblique  shock  without  further  study.  Using  this  first  approx¬ 
imation  as  an  indication,  we  note  also  that  the  maximum  gradient  of 
the  total  density  (emitted,  scattered,  and  incident  families)  grows 
in  an  exponential  fashion  from  the  leading  edge  and  does  not  lie  on 
the  line  x  =  y.  This  is  the  line  which  would  show  up,  for  example, 
on  a  Schlieren  photograph. 
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The  tests  of  Ref.  1  consist  of  2  sets  of  4  runs  each  at  varying 
Mach  numbers,  each  set  at  constant  wall- to- stagnation  temperature. 

From  Eq.  (33)  we  find  that  there  is  to  be  expected  an  asymptotic 

,  Ap /p  * 

limit  of  the  ratio  — —  as  Mq  goes  to  infinity 

o  _ 

T  ,  AP/P  CD  .  F  .  1 7-1  Tw  ,  F./y-l  Tw  _  .... 

Lim  M  2  1  “  C  V  2  T  2  *  2  T  L 

Lim  o  s  J  s 

M  0 
o 


It  is  convenient  to  use  L  in  the  iterative  process,  although  there 
are  no  data  available  to  determine  it  directly  from  the  experiment 
(which  would  be  convenient) . 


We  have 


2L  -  F. 


7-1  Tw 


Ea/7-1  _w 
C  ’  2  T 


Having  two  independent  tests  allows  C  *  D  to  be  eliminated.  It 


follows  that  letting 


7-1  Tw 


(Ti  h 


T2  L2/ 


F  ti  -  F2  +  2 5  -  V 


As  we  stated  in  the  text,  we  begin  by  assuming 


F  =  1;  7  =  1.4 


L  is  used  here  as  a  shorthand  notation,  as  per  Eq.  (A-l) 
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(This  third  constant  could  not  be  unambiguously  found  from  only  two 
test  series.)  There  is  some  question  as  to  whether  7  =  1.4  for 
both  tests,  but  any  departure  could  be  absorbed  in  the  unknown  con¬ 
stant  F.  The  procedure  used  is  this: 

a.  Assume  C  and  (for  the  hot  tests);  solve  for  CD  from 
Eq.  (A-2);  solve  Eq.  (A-l)  at  two  Mach  numbers  and  compare 
with  experiment. 

b.  Solve  for  from  Eq.  (A-3)  (f°r  the  cold  tests);  solve  Eq.  (A-l) 
at  two  Mach  numbers  and  compare. 

Several  tries  lead  to  the  curves  shown  in  Fig.  6,  with  the  values 
F  -  1 
7  «  1.4 
C  -  0.33 
D  *  25 

Effect  of  Changing  F 

Note  that  this  method  of  simultaneous  solution  (Eqs .  (A-l), 

(A-2),  and  (A-3))  makes  the  variables  C,  D,  and  F  interdependent. 

To  demonstrate  this.  Fig.  A-l  shows  the  effect  of  varying  F  at  con¬ 
stant  C.  One  has 

C  -  0.33  -  F  =  1  ,  D  =  25 

C  -  0.33  -  F  =  1.2  D  =  27 

This  new  value  of  F  corresponds  to  thermal  accommodation  less  than 
unity  (approximately  1/F  =  80  per  cent) .  It  is  seen  that  this  improves 
the  simultaneous  correlation  of  theory  and  experiment. 

The  small  change  in  F  (decrease  of  the  accommodation  coefficient) 


did  not  destroy  the  fit  for  the  hot  tests.  However,  were  F  to  be 
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increased  further.,  C  would  need  to  be  decreased.  Numerically,  the 
fit  would  be  improved  for  both  tests,  but  the  term 

I  /y-i Tw 

C  V  2  Ts 


would  keep  increasing  to  values  close  to  unity.  This  invalidates 

the  fundamental  assumption  V  /V  <<:  1  and,  therefore,  the  whole 

e  n  ’  ’ 

analytical  result. 


Discussion  of  Heat-Transfer  Results 

We  turn  now  to  the  heat-transfer  measurements  available  to  date. 
There  are  two  tests  which  show  clearly  the  existence  of  a  heat- 
transfer  plateau  at  about  the  same  distance  from  the  leading  edge 
as  the  pressure  plateau.  Their  magnitudes  are 

Cu  M3  er  140  at  M  =  19.2  T  /T  »  0.083 
H  w  s 

C„  M3  S'  320  at  M  =  25.1  T  /T  -  0.212 
H  w  s 


where  C  is  the  Stanton  number  based  on  the  net  heat  transfer  to  the 
ri 

wall . 


As  outlined  in  the  text  (see  Fig.  8),  we  form  the  ratio  C „/C 

H  HFMl 

with  the  over-all  incident  free-stream  (translational)  energy- transfer 
formula 


0.288 

M 

o 


Personal  communications  from  H.  T.  Nagamatsu. 
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for  y  -  1.4.  This  yields 


"'HFMiy  max 


=  1.24  at  M  =19.2 
o 


c  ,  =  1.76  at  M  =  25.1 

^  HFMi/ nax 


The  theory  would  predict 


where  n  /n  is  a  weak  function  of  the  value  of  C  and  M  , 
s  e  o 

Equation  (A-4)  provides  one  more  relationship  which,  when 
solved  simultaneously  with  Eqs.  (A-l)  through  (A-3),  can  serve  to 
define  unambiguously  the  set  of  parameters  C,  D,  F.  However,  it 
can  be  seen  immediately  from  the  low  value  of  C  /C  .  and  from 
Fig.  8  that  the  solution  of  Eq.  (A-4)  will  yield  value  of  F  too 
large  to  remain  in  the  domain  of  validity  of  this  cold-wall  theory. 

Concluding  Remarks 

In  summary,  the  theoretical  trends  are  found  to  be  in  accord 
with  all  the  trends  exhibited  by  the  experimental  data,  both  those 
relating  to  the  induced  pressure  and  those  concerning  the  leading- 
edge  heat  transfer. 

However,  the  comparison  of  data  with  theory,  as  well  as  the 
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inspection  of  the  heat-transfer  measurements,  indicates  that,  this 
series  of  tests  involves  a  thermal-accommodation  coefficient  which 
is  very  low.  Rough  calculations  show  that  convergence  of  all  tests 
requires  F  of  about  2.5;  that  is,  a  the rma 1 - a c commo da t ion  co e f f i cient 
of  about  40  per  cent. 

This  is  also  an  indication  of  strongly  specular  reflection. 

Consequently,  the  assumption  of  9  =  rr/2  made  in  the  analysis 
* 

should  be  revised. 

Apparently,  in  spite  of  the  low  wall  temperature  the  wall- 
emission  velocity,  Vg,  in  the  experiments  of  Ref.  1  was  not  sufficiently 
low  to  satisfy  the  conditions  implied  in  the  theory.  Although  it 
is  possible  to  find  a  consistent  set  of  empirical  correlation  para¬ 
meters  to  fit  the  data  and  theory  together,  this  is  not  meaningful 
until  more  tests  properly  planned  as  indicated  by  the  theory  become 
available . 


Note  that  the  effect  of  8  can  be  partly  compensated  for  by  choosing 
D  properly.  This  is  the  reason  that  a  better  empirical  correlation 
of  the  pressure  data  with  F  =  1  can  be  obtained  than  of  the  heat- 
transfer  data. 
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